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Abstract 

We present a new heuristic point-to-point routing algorithm based on contraction hierarchies (CH). 
Given an e > 0, we can prove that the length of the path computed by our algorithm is at most (1 + e) 
times the length of the optimal (shortest) path. CH is based on node contraction: removing nodes from 
a network and adding shortcut edges to preserve shortest path distances. Our algorithm tries to avoid 
shortcuts even when a replacement path is e times longer. 

1 Introduction 

Routing in static road networks is essentially solved. There exist fast routing algorithms that are exact. 
However, for other graph classes, e.g. unit disk graphs as they appear in sensor networks, the algorithms 
do not work very well. Also when several objective functions should be supported within a road network, 
current algorithms [4] face some problems since the hierarchy of the network is different for different objective 
functions, e. g. time and distance. One possibility to alleviate these problems is to drop the exactness of the 
algorithms and allow some error. We show how to adapt contraction hierarchies (CH) [5] so that we can 
guarantee an error of e. 

Related Work 

There has been extensive work on speed-up techniques for road networks. We refer to [2, 6, 1] for an overview. 

We can classify current algorithms into three categories: hierarchical algorithms, goal-directed approaches 
and combinations of both. Our algorithm is based on contraction hierarchies (CH) [5] . A CH uses a single 
node order and contracts the nodes in this order. A slightly modified bidirectional Dijkstra shortest path 
search then answers a query request, touching only a few hundred nodes. For our algorithm, we modify the 
node contraction, i.e. the decision which shortcuts we have to add, and the query. 

Recently, also speed-up techniques for time- dependent road networks have been developed ([3] presents 
an overview). 

Also algorithms to support multiple objective function have been developed lately. The most recent one 
is an adaption of the CH algorithm [4] to allow a flexible query. 

2 Heuristic Node Contraction 

CH performs precomputation on a directed graph G = (V, E) with edge weight function c : E K. + . The 
precomputation of CH is based on node contraction: a node u is contracted by removing it from the network 
in such a way that shortest paths in the remaining graph are preserved. In the exact scenario, where we 
want to find the optimal path, we use Algorithm 1. To preserve all shortest path distances, it is enough to 
preserve the shortest paths distances between the neighbors of u. So given two neighbors v and w with edges 
(v, it) and (u, w), we should find the shortest path P between v and w avoiding u (Line 4). When the length 
of P is longer than the length of the path (v,u, w), a shortcut between v and w is necessary with weight 
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c(v, u) + c(u, w). When the length of P is not longer, no shortcut is necessary and P witnesses this, in this 
case we call P a witness. 

In the heuristic scenario, we use Algorithm 2. Intuitively, we want to avoid more shortcuts than in the 
exact scenario, i.e. a shortcut between v and w should be avoided, even when the path P is a bit longer 
than (v, u, w). To be able to guarantee a maximum relative error of e, we need to ensure that the errors 
do not stack when a node on P is contracted later. So we store another edge weight c with each edge. 
Intuitively, c is the witness memory of an edge. When a witness P prevents a shortcut, even though in the 
exact scenario, the shortcut would be necessary, the witness must remember this delta. We let the edges of 
the witness remember this (Lines 9-11), so that c(P) < c(v,u) + c(u,w). Note that 7 < e (Line 9), which 
gives Lemma 1. 



Lemma 1 For each edge (v,w) holds 



c(v, w) 
1 + e 



< c(v, w) < c(v, w) 



Note that there are other ways to ensure that the errors do not stack. Currently, we proportionally 
distribute the delta among all edges of the witness. We could distribute the delta differently or even try to 
find other potential witnesses. Also, avoiding a shortcut can lead to more shortcuts later, as every shortcut 
is a potential witness later. 

Algorithm 1: SimplifiedExactConstructionProcedure(G = (V,E),<) 



foreach u G V ordered by < ascending do 
foreach (v, u) G E with v > u do 
foreach (u, w) G E with w > u do 

find shortest path P = (v, . . . , w) $ u; 
if c(P) > c(v, u) + c(u, w) then 
[_ E:= E U {(v, w)} (use weight c(v, w): 



c(v, u) + c(u, w)) 



Algorithm 2: SimplifiedHeuristicConstructionProcedure(G = (V, E),<,s) 



1 c:= c; 

2 foreach u £ V ordered by < ascending do 



// store second weight per edge 
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foreach (v, u) 6 E with v > u do 
foreach (u, w) G E with w > u do 

find shortest path P = (v, . . . , w) ^ u; 
if c(P) > (1 + s)(c(v, u) + c(u, w)) then 

I E:= E U {(v, w)} (use weight c(v, w):= c(v, u) + c(u, w), c(v, w):= c(v, u) + c(u, w)); 



else 



' ' c(v,u)-\-c(u,w) ' 



// c(P) = (1 + j)(c(v, u) + c(u, w)) 



foreach (x,y) G P do 
c(x,y):= min |c(x,y), 



3 Heuristic Query 

The basic query algorithm is the same as for exact CH (Algorithm 3). It is a symmetric Dijkstra-like 
bidirectional procedure. We will introduce it formally and prove its correctness. 

The query algorithm does not relax edges leading to nodes lower than the current node. This property is 
reflected in the upward graph G^:— (V,E^) with E^:= {(u,v) G E \ u < v} and, analogously, the downward 
graph G 4 := (V,^) with E±:= {(«,«) G E \ u > «}). 



We perform a forward search in G^ and a backward search in G^. Forward and backward search are 
interleaved, we keep track of a tentative shortest-path length and abort the forward/backward search process 
when all keys in the respective priority queue are greater than the tentative shortest-path length (abort-on- 
success criterion). 

Both search graphs G-f and Gj. can be represented in a single, space-efficient data structure: an adjacency 
array. Each node has its own edge group of incident edges. Since we perform a forward search in G-|- and a 
backward search in Gj., we only need to store an edge in the edge group of the lower incident node. This 
formally results in a search graph G* = (V, E*) with £j := {(v, u) | (u, v) G E±} and E* := EfUE^. Finally, 
we introduce a forward and a backward flag such that for any edge e G E* , t (e) = true iff e G E-f and 
4- (e) = true iff e G E^. Note that G* is a directed acyclic graph. 

Algorithm 3: ExactQucry(s,t) 

1 dy.= (oo, . . . , oo); d-f [s]:= 0; d(.:= (oo, . . . , oo); dj.[i]:= 0, d:= oo; // tentative distances 

2 Q t = {(0,s)}; Qi = {(0,t)}; r:= f; // priority queues 

3 while (Q-j- ^ or Qi ^ 0) and (d > min {min Qf, minQ^}) do 

4 if Q~, r ^ then r:= ^r; // interleave direction, -i J [=\. and -i L=t 

5 (-,w):= Q r .deleteMin(); d:= min {d, d^[u] + tZJu]}; // u is settled and new candidate 

6 foreach e = (u,v) £ £* do // relax edges of u 

7 if r(e) and (d,[u] + c(e) < cf r M) then // shorter path found 

8 d r [v]:— d r [u] + c(e); II update tentative distance 

9 Q r .update(d r // update priority queue 

10 return d; 



Lemma 2 Let P be an s-t-path P in the heuristic CH (G = (V,E), <,e). Then there exists an s-t-path P' 
in the CH of form (s = uq, u\, . . . , u p , . . . , u q — t) with p, q £ N, Ui < Uj+i for i G N, i < p and Uj > %+i 
for j £N,j) < j < q, short path form (PF) And for P' holds c(P') < c(P). 

Proof. Given a shortest s-t-path P = (s = uq, Ui, . . . , u p , . . . , u q — t) with p, q G N, u p — maxP, that is 
not of the form (PF). Then there exists a k G N, k < q with Uk-i > Ufe < Uk+i- We will recursively construct 
a path of the form (PF). 

Let Mp:= {uk \ Ufe-i > Uk < u^+i} denote the set of local minima excluding nodes s, t. We show that 
there exists a s-t-path P' with Mp* = or minA/p/ > min Mp. 

Let Mfe:= min Mp and consider the two edges (uk-i,Uk), {uk,Uk+i) G P. Both edges already exist at the 
beginning of the contraction of node Uk- So there is either a witness path Q — (uk—i, ■ • ■ , Ufc+i) consisting of 
nodes higher than uu with c{Q) < (1 + e)(c(Uk-i, Wfe) + c{uu, Wfe+i)) or a shortcut (Ufc-i, u/s+i) of the same 
weight is added. So the subpath -P| Ufc _i-s-« fe+ i can either be replaced by Q or by the shortcut (uk-i, Uk+i)- 
If we replace the subpath by Q, our construction ensures that c(Q) < c(P| Ufc _ 1 _> Ufe+1 ). Also if we added 
a shortcut, c{uk—x, Uk+i) < 2(P| Mfc _ 1 ^ Ufc+1 ) holds. So the resulting path P' consists of nodes higher than 
Uk and has the property c(P') < c(P). Since n < oo, there must exist a s-t-path P" with Mp>/ = 0, 
c(P") < c(P) and of the form described in (PF). □ 

Theorem 1 proofs the correctness of our basic query algorithm. 

Theorem 1 Given a source node s and a target node t. Let d(s,t) be the distance computed by the heuristic 
CH algorithm with e > and let d(s,t) be the optimal (shortest) distance in the original graph. Then 
d{s,t) < d{s,t) < (1 + e)d(s,t). 

Proof. Let (G = (V,E),<,e) be a heuristic CH with e > 0. Let s,t G V be the source/target pair 
of a query. It follows from the definition of a shortcut, that the shortest path distance between s and t 
in the CH is the same as in the original graph. So d(s,t) < d(s,t) holds. Every shortest s-t-path in the 
original graph still exists in the CH but there may be additional s-t-paths. However since we use a modified 
Dijkstra algorithm that does not relax all incident edges of a settled node, our query algorithm does only 



find particular ones. In detail, exactly the shortest paths of the form (PF) are found by our query algorithm. 
From Lemma 2, we know that if there exists a shortest s-t-path P then there also exists an s-i-path P' of 
the form (PF) with c(P') < c(P). Because of Lemma 1, we know that -fTF — c(-P') and c(P) < c(P) so 
that c(P') < (1 + e)c(P). So our query algorithm will either find P' or another path, that is not longer than 
P' . □ 
Although we use c in the proof, the query algorithm does not use it at all. So we only need to store it 
during precomputation but we do not need to store it for the query. 

4 Heuristic Stall-on-Demand 

In the previous section, we proved that the basic heuristic query algorithm does not need any changes 
compared to the exact scenario. However, an important ingredient of a practically efficient query algorithm is 
the stall-on-demand technique first described by Schultes et al. [7] . This single improvement brings additional 
speedup of factor two or more. 

While performing a query, we do not relax all edges of a node, we only relax edges leading upwards. Here 
and in the following we well assume that the query is from source node s to target node t and we will only 
consider the forward search. The backward search is completely symmetric to the forward search. We call 
a path leading only upwards being an upward path. In the exact scenario, we check, whether a settled node 
u in the forward search is reached suboptimally via a s-u-path P = (s = Vi, . . . , = u), i.e if there is a 
s-u-path P' with c(P) > c(P'). If this was the case, we can safely stop (stall) the search at u since it would 
never contribute to a shortest path. The path P is an upward path and stalling happens if and only if the 
path P' is not an upward path.. 

However, in the heuristic scenario with e > 0, we would destroy the correctness of our algorithm when 
we would apply the same rule. Consider as example the graph in Figure 1. During the contraction of u, no 
shortcut for the path (x, u, v) is added since the path (x, y, v) is a witness that is just a factor 1 + s larger. 
The forward query starting at s should settle the nodes in the order s, x, y, u, v, z. However, if we would not 
change the stalling condition, we would stall u while settling it because the path (s,u) is longer than the 
path (s, x, u), which is not an upward path. Furthermore, we would propagate the stalling information to v, 
so node v reached via upward path P — (s, x, y, v) gets stalled by the shorter path P' = (s, x, u, v). So node 
v is stalled and we would never reach node z with the forward search and therefore could never meet with 
the backward search there. 




Figure 1: The stalling condition of the exact query fails in a heuristic CH, as node z is never reached in the 
forward search from s. 

To ensure the correctness, we change the stalling condition. We split the path s-u-path P' in a path P[ 
and P2 so that P[ is the maximal upward subpath. Let x be the node that splits P' in these two parts. 
Then we stall u only if node x is reached by the forward search and 



c(P{) + (l + s)c(Pi)<c(P) 



(1) 



The symmetric condition applies to the backward search. 



Algorithm 4: ExactQucry(s,i) 



1 df.= (oo, . . . , oo); d^[s):— 0; <i^:= (oo, . 

2 Q t = {(0,s)}; Q ; = {(0,t)}; r:= f; 



. ,oo); d4.fi] := 0, d:= 00; 



3 while (Q t ^ or Q4. ^ 



if Q^ r ^= then r:= -^r; 
(-,u):= Q r .deleteMin(); d:= min {d, df [u] - 
if isStalled(r, u) then continue; 
foreach e = (u, v) £ P* do 

if r(e) and (d r [u] + c(e) < c? r [f]) then 
d r [«]:= c? r [u] + c(e); 
Qr-update((i r [i»],u); 
if isStalledfr, v) then unstall(r, w): 



// tentative distances 
// priority queues 
and (d > min {min Q^, min Qi}) do 

// interleave direction, -1 t=i an d ^ 4-=t 
g^[u]}; // 11 is settled and new candidate 

// do not relax edges of a stalled node 
// relax edges of u 
II shorter path found 
// update tentative distance 
// update priority queue 



if (-ir)(e) and (d r [v] + c(e) < then // path via v is shorter 

stall(r, u, d r [v] + c(e)); // stall u with stalling distance d r [v] + c(e) 

break; // stop relaxing edges of stalled node u 



15 return d; 



Algorithm 5: HcuristicQuery(s,i) 



1 d^:— (00, . . . , 00); d-j-[s]:= 0; d±:= (00, . . . , 00); dj,[i]:= 0, d:= 00; 

2 Q t = {(0,s)}; Q^ = {(0,i)}; r:= t; 

3 while (Q-j- 7^ or Q4 7^ 0) and (d > min {min Q-)-, min Qi}) do 



// tentative distances 
// priority queues 
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if Q_, r ^ then r:= ^r; // interleave direction, -1 t=| and -1 J.=j" 

(•,it):= Q r .deleteMin(); d:= min {d, d-j-[u] + d^[u]}; II u is settled and new candidate 

if isStalled(r , u) then continue; // do not relax edges of a stalled node 

foreach e = (u,v) £ E* do // relax edges of u 

if r(e) and (d r [u] + c(e) < d r [u]) then // shorter path found 

d r [v]:= d r [u] + c(e); II update tentative distance 

Q r .xvpda,te(d r [v],v)] II update priority queue 
if isStalledfr, v) then unstall(r, v); 

if (-ir)(e) A d r [v] + (1 + s)c(e) < d r [u] then // path via v is shorter 

stall(r, u, d r [v] + (1 + e)c(e)); // stall w with stalling distance d r [v] + (1 + e)c(e) 
break; // stop relaxing edges of stalled node u 



15 return d; 



To prove that the stall-on-demand with (1) is correct, we will iteratively construct in Lemma 3 a new 
path from a stalled one. 

Lemma 3 Let (P,v,w) be a stall state triple (SST): P being an s-t-path of the form (PF), node v being 
reached by the forward search by P\ s ^ v o-nd not stalled and node w being reached by the backward search by 
P\ w^rt and not stalled. Define a function g on an SST: 

g(P, v, w) := c(P| s _h,) + (1 + e)c(P| ) + c(P\ 

If one of the nodes in P\ v ^ w becomes stalled, then there exists an SST (Q 7 x,y) with 

g(Q,x,y) < g{P,v,w). 

Proof. Let u £ P\ v ^ w be the node that becomes stalled. WLOG we assume that P\ v ^ u is an upward 
path, i.e. the stalling happens during the forward search. Then there exists an s-w-path P' that is split in P[ 



and P' 2 as defined in (1) so that c(P{) + (1 + e)c(P2) < c(P\ s ^ u ). Let x be the node that splits P' into these 
two subpaths. Let R be the path of form (PF) that is constructed following Lemma 2 from the concatenation 
of P' 2 and P\ u ^ w . Let Q be the concatenation of P[, R and P\ w ^t and y:= w. By construction, (Q, x, y) is 
a SST and we will prove that it is the one that we are looking for. 

g(Q,x,y) = c(Q\ s ^ x ) + (1 + e)c(Q\ x ^ y ) + c(Q\ y ^ t ) 

d ^ c(P{) + (1 + e)c(R) + c(P\ u ^ t ) 

< 2 c{P{) + (1 + e)(c(i3) + c(P|„^,)) + c(PU) 

L < c{P[) + (1 + £ )c(P 2 ') + (1 + e)c{P\ u ^ w ) + c(P|^ t ) 
(i) 

< c(P| s ^„) + (1 + e)c{P\ u ^ w ) + c{P\ w ^ t ) 

= c(P\ s ^ v ) + c(P\ v ^ u ) + (1 + e)c(P\ u ^ w ) + c(P\ w ^ t ) 

L < c(P| s ^) + (1 + e)c(P\ v ^ u ) + (1 + £ )c(P| u ^) + c{P\ w ^ t ) 
= 9{P,v,w) 

□ 

With Lemma 3 we are able to prove the correctness of heuristic stall-on-demand (1) in Theorem 2. 

Theorem 2 Theorem 1 still holds when we use heuristic stall- on- demand (1). 

Proof. The proof will iteratively construct SST's with Lemma 3 starting with the path P found in the 
proof of Theorem 1/Lemma 2 and the nodes s and t. Obviously, at the beginning of the query, both nodes 
s and t are reached and not stalled, so (P, s, i) is an SST and 

g(P, s, t) = c(P| s ^ s ) + (1 + e)c(P|._, t ) + c(P| t _rt) = (1 + e)c(P) < (1 + e)d(s, t)- 

We will prove that after a finite number of applications of Lemma 3, we obtain an SST (Q,x,y) so that Q 
is found by our query with stalling. For this path Q holds: 
c{Q) = c(Q| s _>. x ) +c(Q\ x ->y) +c(Q\ y -n) 

L.l 

< c{Q\ s ^ x ) + (1 + e)c(Q\ x ^ y ) + c(Q\ y ^t) 
= g{Q,x,y) 

L.3 

< g{P,s,t) 

< (l + e)d(s,t) 

Since our graph is finite, and due to the "<" in Lemma 3, we can apply Lemma 3 only finitely many 
times. The final SST (Q,x,y) will be found by our query algorithm since x is is reached in the forward 
search and not stalled, y is reached in the backward search and not stalled. And since this is the final SST, 
no node on the path Q\ x ^ y will be stalled so our query will find the path Q or a shorter path. □ 



5 Applications 

We described our heuristic for just a single edge weight function. Further applications may be time-dependent 
road networks, where not only the travel time functions are approximations, but also the shortcuts. Also, 
it can help for multi-criteria optimization. It would be simple to extend it to the flexible scenario [4] with 
two edge weight functions. This would hopefully make the splitting of parameter intervals unnecessary, as 
approximation can reduce the number of shortcuts. Without splitting, preprocessing time and space are 
greatly reduced, this comes in handy as they are significantly larger than in the single-criteria scenario. 

In general, heuristic CH can help to process graphs that are less hierarchically structured than road 
networks. In theory, every graph can be preprocessed by CH. But in practice, this is not feasible. The 
witness searches that decide the necessity of shortcuts take too long, especially when the remaining graph 
is very dense after a lot of shortcuts have been added. Reducing the number of shortcuts can therefore 
significantly reduce the preprocessing time. 



6 Conclusion 



We developed an algorithmic idea for an heuristic CH algorithm that guarantees the approximation ratio. 
Now experiments have to show the feasibility of this approach in practice. 
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